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Abstract 


In the present work a theoretical model of a solid oxide electrolyzer based on an electrolyte having both oxygen ion and proton conductivity is 
considered. The main parameters of the electrolytic process and an electrolyzer (distribution of gas components, electromotive forces and current 
densities along the electrolyzer channel, average values of electromotive forces and current densities) were calculated depending on a proton 
transport number and mode of the reactants’ feeding (co- and counter-flow). The performed analysis demonstrates considerable influence of the 
mode of feeding on all parameters of the electrolyzer: operation under the counter-flow mode is preferable as regards the specific characteristics 
and uniformity of their distribution along the electrolyzer. It is shown that the electrolyser’s specific characteristics increase with the increase of 


the proton transport number. 
© 2007 Elsevier B.V. All rights reserved. 


Keywords: Solid oxide electrolyzer; Co-ionic electrolyte; High temperature steam electrolysis 


1. Introduction 


Hydrogen production by means of water electrolysis is a very 
attractive way for pure hydrogen production by using electri- 
cal energy during a period when using the electricity by other 
consumers is minimal. Advanced low temperature electrolyz- 
ers demonstrate rather low power inputs, about 4.1 kWh/m? Ho, 
which corresponds to an average cell voltage of 1.71 V. 

An idea to use electrochemical devices based on solid oxides 
for hydrogen production by means of steam decomposition 
appeared in the late 1960s [1]. In the late 1970s Rohr [2] per- 
formed series of experiments with electrolysis cells and stacks 
and showed that they have acceptable stability during ca. 8000 h. 
Since then the base of the electrolysis cells was yttria-stabilized 
zirconia (YSZ) having oxygen ion conductivity. It was demon- 
strated that a high temperature solid oxide electrolyzer (SOE) 
based on an oxygen ion solid oxide electrolyte had acceptable 
characteristics when it worked at 1.1 V [3], which corresponded 
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to electrical power inputs 2.63 kWh/m? Hp. In this case it is 
necessary to supply the SOE electrochemical section and the 
vaporizer with additional heat (about 0.5 kWh/m? H; each), but 
both electrical power inputs and total (electrical + heat) inputs 
are significantly lower than in the low temperature electrolyzer. 

In the beginning of 1980s, oxides having proton conductivity 
were discovered and at the very early stage it was proposed to 
use them for steam electrolysis [4,5]. Since then a number of 
experiments were performed with the cells based on protonic 
electrolytes (see, for instance, [6,7]). Recently, the interest to 
the application of protonic electrolytes for steam electrolyzers 
is growing [8,9]. 

Itis known now that only a little part of “protonic electrolytes” 
has practically pure proton conductivity. A great number of these 
solid oxides possess a significant part of oxygen ion conductiv- 
ity in addition to proton conductivity, so it is logical to call them 
“co-ionic electrolyte” [10]. It was shown that an achievable effi- 
ciency of a solid oxide fuel cell (SOFC) increases, as the proton 
transport number of the electrolyte increases [11]. However the 
most known proton electrolytes are unstable in CO2-containing 
atmosphere which is typical to the SOFC. The problem of car- 
bonates formation is absent in the SOE for steam splitting, so it is 
possible to use a co-ionic electrolyte in the SOE. It is interesting 
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Nomenclature 

a anode 

c cathode 

E electromotive force 

F Faraday constant 

h electrolyte thickness 

H enthalpy 

i current density 

I current 

j mass flux 

K equilibrium constant of the reaction of hydrogen 
oxidation 

L length of the channel-working zone 

p dimensionless partial pressure; heat flux 

R gas constant 

t transport number 

T absolute temperature 

U voltage 

x coordinate of the channel working zone 

X component of the gas mixture 

Greek letter 

o conductivity 

Subscripts 

eff effective value 

F related to Faraday process 

H related to proton transport 

H2 hydrogen 

H20 steam 

(0) related to oxygen ion transport 

O2 oxygen 

th thermoneutral 

tot total 

T absolute temperature 

v related to convective mass transport 


from theoretical point of view to analyze the influence of ion 
transport numbers on the SOE characteristics. 


2. Theoretical model 


A simplified scheme of the processes occurring in the SOE 
based on the co-ionic electrolyte is presented in Fig. la. The SOE 
consists of the electrolyte membrane opposite sides of which 
are in contact with porous electrodes, an anode and a cathode, 
which, in turn, are in contact with gas flows moving within corre- 
sponding gas channels. The electrodes are connected with a DC 
power source. In general case, steam is fed to both the cathode 
and anode channel. 

At the anode, water is decomposed to hydrogen and oxygen 
ions according the following electrochemical reaction: 


H20 + 2e7 = H +0% (1) 


Fig. 1. A schematic diagram of the processes in the SOE based on aco-ionic elec- 
trolyte (a); an electrolytic cell based on a co-ionic electrolyte (b); an equivalent 
circuit of the cell (c). 


At the cathode, decomposition of water occurs as follows: 
H20 = 0.502 +2H* + 2e7 (2) 


Protons move across the electrolyte to the anode where they 
form hydrogen: 


2H* + 2e7 = H, (3) 


So, reactions (1) and (3) are two ways for hydrogen formation 
at the anode of the SOE. Oxygen ions forming at the cathode 
according to the reaction (1) move to the cathode where they 
form oxygen: 


O% = 0.502 + 2e7 (4) 


This reaction and the reaction (2) are mechanisms of oxygen 
formation in the SOE. 

From the electrochemical point of view, the SOE presents an 
ionic conductor opposite sides of which are in contact with dif- 
ferent gas atmospheres (Fig. 1b). The electrolyte has two types 
of charge carriers, namely, oxygen ions and protons and can 
be characterized by corresponding conductivities, co and oy. 
Each gas phase can be characterized by the dimensionless partial 
pressures of oxygen, hydrogen and steam, po,, pH, and ppo, 
respectively. The first two values are equal to the corresponding 
mole fractions of the components if they are mixed with steam; 
in case they are chemically bonded, the corresponding partial 
pressures are calculated by using the following equations: 


2 
PH20 
o = | —— (5) 
i (ge) 
PH20 
pm = 23S, © 
Kpo, 


where K is the equilibrium constant of the reaction of hydro- 
gen oxidation. 

An equivalent circuit of the SOE based on the co-ionic elec- 
trolyte presents two parallel galvanic cells connected with an 
external current source (Fig. 1c). Each galvanic cell is due to 
presence of the given charge carrier and is characterized by an 
electromotive force, EMF, due to drop of corresponding partial 
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pressure, and by an internal resistance to transfer of correspond- 
ing charge carrier. The electromotive forces are calculated by 
the following equations: 


RT AE- 
Eo =—-— 1n POs (7) 
4F Po 
RT A 
ea mne (8) 
2F Pup 


where R, T, F are the ideal gas constant, absolute temper- 
ature and Faraday’s constant, respectively. The resistances are 
calculated as follows: 


R 

Ro = — (9) 
to 
R 

Ry = — (10) 
tH 


where fo, ty are transport numbers of oxygen ions and pro- 
tons, respectively, determined as follows: 


90 


to = (11) 
oO 

co (12) 
oO 


In the frame of this consideration, it is assumed that resis- 
tances of the electrochemical reactions (1)—(4) are negligible as 
compared with the resistances for ions transfer (9) and (10). 

Applying laws of electrical circuits to the circuit of Fig. Ic 
partial currents connected with transport of corresponding ions, 
Io and Fy can be calculated. It is more convenient to use values 
of current densities, ig and iy, which can be presented by the 
following equations: 


ne moo) (13) 
= am) (14) 


where Ais the electrolyte thickness. Combination of Eqs. (13) 
and (14) gives the total ion current density, itot- 


o(U — E 
(ee (15) 
h 
where 
Eerp = to Eo + tu Eu (16) 


is an effective electromotive force of the cell based on the 
co-ionic electrolyte. 

Fluxes of the components from (to) the electrodes due to 
reactions (1)—(4), called Faraday fluxes, jp, are equal 


i io + iy 

+= 17 
JFM oF (17) 
FH ,0(C) = 2 (18) 
JF,H2,0(C) = OF 


at the cathode and: 


io + in 

= 19 
JF.O> ZF (19) 
l iH 
jFmola) = E (20) 


at the anode; minus means disappearance of the component. 

Gas mixtures moving along the channels can be charac- 
terized by convective fluxes of the components, jy,0,, jv,Ho; 
Jv.Ho0(C), jv,m o(a), and by total fluxes, jy(c) and j,(a), for the 
cathode and for the anode channels, respectively. It is obvious 
that the ratio of the convective flux value of a given component 
to the total convective flux value for a given cross-section of 
the channel just equal to the partial pressure of the component 
in this cross-section. Change of the convective fluxes along the 
channels due to the Faraday fluxes obeys to equations: 


djy H2 M 
H 21 
ay JFM (21) 
djy O2 r 
%2 22 
dz JF, (22) 
djy, m, o(c) . 
O = jru0©) (23) 
dx 
djuwo@ . 
OS, et) (24) 
dx 


The following SOE characteristics were calculated: relative 
fluxes of steam at the inlet of the cathode and anode channels; 
distributions of partial pressures of the components, EMF values 
and partial current densities along the SOE; average values of 
EMFs and partial current densities within the SOE. Each set of 
data was calculated for a number of the proton transport number. 

The following algorithm was used for the calculations. The 
length of the working part of the SOE was set as 1. It was 
divided into great number of intervals (up to 2000). Chosen val- 
ues of partial pressures of the components in the input and output 
gases were assigned. Initial values of inlet convective fluxes to 
the anode and to the cathode channels were assigned as a first 
approximation. From the parameters of the inlet gas mixtures, 
values of EMFs and partial current densities in the initial point 
were calculated using Eqs. (7) and (8) and Eqs. (13) and (14), 
respectively. By using Eqs. (13) and (14) and Eqs. (17)-(24), 
changes of convective fluxes of each component X within the first 
interval were calculated according to the following equation: 


Aju,x = JF x ^x (25) 


Values of the components’ convective fluxes at the end of 
the first interval were calculated and then EMFs and partial cur- 
rent densities in this point were also calculated. Afterwards, the 
procedure was repeated for the next interval, and again right to 
the end of the channel. Obtained partial pressures of the compo- 
nents in the outlet gases were compared with the chosen assigned 
values; in case they differ from each other the initial values of 
the inlet convective fluxes were corrected and the calculations 
were repeated from the beginning. The cycles of the calculations 
were repeated until close agreement between the obtained and 
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the assigned partial pressures of the component was achieved. 
The calculations were performed in the shell Delphi 6.0. 

In order to reveal the main peculiarities of the SOE based on 
a co-ionic electrolyte, the following parameters were assigned. 
The inlet cathode gas: H20+H2, py,o = 0.98; the inlet 
anode gas: H20 +02, pmo = 0.995; the outlet cathode gas: 
H2 + H20, puro = 0.1; the outlet anode gas: O2 + H20, pu,o = 
0.1; T=1100K; modes of gas feeding: co-flow and counter- 
flow; applied voltage: 1.1 V; electrolyte specific conductivity: 
o/h=1Smcm~?. The latter parameter is real for supporting 
electrolyte; its value does not influence the distribution of gas 
components’ concentrations and EMFs along the channels. The 
presented current densities are called “conditional current densi- 
ties”, i.e., corresponding to the condition o/h= 1 Smcm~?. The 
obtained results for the current densities can be easily recal- 
culated for other values o/h. The range of the proton transport 
number was chosen from 0.02 to 0.98. 


3. Results and discussion 


It is clear that the input steam flows to the cathode and anode 
channels depend on the proton transport number: the higher ty 
value the less the steam flow to the cathode channel; on the con- 
trary, the anode steam flow increases with the ty value increase. 
It was found that for a given ty value the inlet flows both to 
the cathode and anode channel practically do not depend on the 
mode of reagents’ feeding. The cathode input flow decreases 
practically linearly with the fy value increase. It is illustrated 
in Fig. 2 where the input flows are presented in relative values: 
they are referred to the unit of produced hydrogen. The anode 
input flow demonstrates also practically linear increasing with 
the ty value increase. 

The results on distributions of the cell characteristics are pre- 
sented for dimensionless relative cell length, x/L, where x is the 
current coordinate of the cell, L is the cell length. It is assumed 
that the point x = 0 corresponds to the cathode inlet. 

Figs. 3 and 4 illustrate the distributions of the partial pressures 
of the electrolysis products, hydrogen and oxygen, along the 
channels for three cases: (i) the electrolyte has practically pure 
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Fig. 2. Dependence of the relative input steam flows on the proton transport 
number. 


oxygen conductivity (t =0.02), (ii) partial conductivities are 
equal (ty =0.5) and (iii) the conductivity is practically protonic 
(ty =0.98). One can see that both the mode of reagents’ supply 
and the proton transport number value influence the distribution. 
The difference of oxygen partial pressures distributions for the 
co- and counter-flow modes is quite obvious first of all due to 
opposite directions of the anode flows. However, it is easy to see 
that the distributions are not symmetrical. 

The proton transport number value also influences strongly 
the distributions of the components in the channels. Rapid 
increase of the oxygen partial pressure at the initial part of the 
anode channel at tg =0.02 is connected with the predominant 
transport of oxygen through the electrolyte into the anode com- 
partment where oxygen is mixed with the small steam flow. On 
the contrary, at ty = 0.98 the steam flow is many times bigger and 
the oxygen partial pressure rises slowly (Fig. 3, left). The differ- 
ence in the rate of the change of the hydrogen partial pressure 
depending on the ty value can be explained analogically. The 
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Fig. 3. Distribution of the oxygen partial pressure along the SOE for the co-flow (left) and counter-flow (right) modes. Number at the lines denotes proton transport 


number. 
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Fig. 4. The same as in Fig. 3 for the hydrogen partial pressure. 
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Fig. 5. Distribution of Ey along the SOE under co-flow (left) and counter-flow (right) mode. Number at the lines denotes proton transport number. 


distributions of the components in the channels allow finding 
EMPs’ distributions. 

Figs. 5 and 6 show distributions of Ey and Eo, respectively, 
for both modes of reagents’ supply. The most remarkable fea- 
tures are several ones. First, under the counter-flow mode, both 
Ey and Eo are lower on average than under the co-flow mode. 
Second, under the counter-flow mode, distributions of both Ey 
and Eo are more uniform than under the co-flow mode. Third, 


1100 


800 


under the co-flow mode, both Ey and Eo change monotonously 
with the fy value increases, the former increases and the latter 
decreases. On the contrary, under the counter-flow mode, both 
Ey and Eo have evident minimum at some intermediate values 
of ty. 

In Figs. 7 and 8 are shown the distributions of the condi- 
tional partial current densities along the SOE for the proton and 
oxygen ion charges, respectively. It can be distinguished that 
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Fig. 6. The same as in Fig. 4 for Eo. 
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Fig. 7. Distribution of the conditional proton current density along the SOE under co-flow (left) and counter-flow (right) mode. Number at the lines denotes proton 


transport number. 


under the counter-flow mode, both the proton and the oxygen 
ion current densities are more uniform and slightly higher than 
under the co-flow mode. The obtained results are of the great 
importance for the calculations of the partial convective flows 
of the reagents and products in both the anode and the cathode 
channels. Besides, they allow calculating the distribution of heat 
flux, po, along the SOE by using equation: 


Po = ito(U — Un,r) (26) 
where Ujp,7 is thermoneutral voltage determined as: 

mas 27 
i 2F 


where AHF is the enthalpy change of the hydrogen oxidation 
reaction at temperature T. In the frame of this work the distri- 
bution of the heat flux was not calculated. However, the main 
conclusion is obvious: the heat flux distribution is more uniform 
under the counter-flow mode; hence temperature distribution in 
the SOE in this case must also be more uniform. 

The advantages of the counter-flow mode of the reagents’ 
supply as well as the use of the co-ionic electrolyte with higher 
ty value become more obvious when analyzing Fig. 9. The most 
important characteristic for the SOE on the whole is the average 
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effective EMF value which can be found by numerical integra- 
tion of the Eo and Ey distributions and using Eq. (16). It is 
interesting that under the co-flow mode, the average effective 
EME has a maximum at fy ~ 0.3 whereas the average effec- 
tive EMF is minimal at ty ~ 1 under the counter-flow mode. 
This means that the region of ty % 0.3 is unfavourable for the 
co-flow mode of operation because in this case the generated 
current density is minimum. For both modes of operation, the 
proton electrolyte demonstrates the minimum average effective 
EMF and, hence, the maximum total current density can be 
achieved. 

Fig. 10 confirms the said above. Presented here are the aver- 
age values of the partial and total current densities obtained 
by numerical integration of the current densities distributions 
depending on the proton transport number. On the whole, the 
average value of the total current density is visibly higher under 
the counter-flow than under the co-flow mode. It is necessary 
to point out that the difference increases with the decrease of 
the applied voltage that can be achieved in thin-film design 
of the SOE. The value o/h for thin-film systems can be sig- 
nificantly higher. For instance, according to [12] the value 
o/h 10S cm’. If such electrolyte will be used in the SOE then 
the average current can reach 0.95 A cm7? under counter-flow 
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Fig. 8. The same as in Fig. 6 for the oxygen ion current density. 
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Fig. 10. Dependence of the average partial and total conditional current densities on the proton transport number under co-flow (left) and counter-flow (right) mode. 


mode and 0.6 A cm7? under co-flow mode under the applied 
voltage U = 1.05 V. 


4. Conclusion 


The performed analysis shows that the co-ionic electrolyte 
can be successfully utilized as a base of solid oxide electrolyzer. 
One of interesting distinctions of this electrolyzer from the one 
based on the oxygen ion electrolyte is the mode of feeding by 
steam: the latter must be fed to both the anode and cathode 
channels, the part of the cathode feeding being higher with ty 
increasing. Another interesting peculiarity of the electrolyzer 
based on the co-ionic electrolyte is more uniform distributions 
of its specific characteristics (current density and specific heat 
flow rate) under the counter-flow mode of feeding. It is interest- 
ing also that the average current density is visibly higher when 
the SOE operates under the counter-flow mode comparing to 
operation under the co-flow mode, the difference being higher 
with the applied voltage decreasing. 

The most promising application of the SOE is its utilization 
for hydrogen production when it is supplied by electrical energy 
in combination with high potential heat for heating the elec- 
trochemical section and with low potential heat for heating the 
vaporizer. The mentioned heat can be supplied from a nuclear 
power station or from a concentrator of solar energy. In this case, 


the considerable part of expensive electrical energy is replaced 
by cheap heat energy; this is impossible for low temperature 
electrolyzers. The SOE based on the co-ionic electrolyte has 
higher specific characteristics comparing to the SOE based on 
the oxygen ion electrolyte. This is a base either for increase of the 
SOE efficiency under the same size of the unit, or for decrease of 
the unit size under the same efficiency. So, utilization of the co- 
ionic electrolyte in the SOE is very prospect and development 
of the co-ionic electrolytes is very actual. 
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